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Abstract 

(— i | Spin measurements are crucial in distinguishing major scenarios of TeV scale new physics 

once it is discovered at the LHC. We give a brief survey of methods of measuring the spin of new 
physics particles at the LHC. We focus on the case in which a long lived massive neutral particle 
is produced at the end of every cascade decay and escape detection. This is the case for Im- 
parity preserving supersymmetry, Little Higgs models with T-parity, extra-dimensional models 
with KK-parity, and a large class of similar models and scenarios. After briefly commenting on 
measuring spin by combining mass and rate information, we concentrate on direct measurement 
■ by observing angular correlations among decay products of the new physics particles. We survey 

[ a wide range of possible channels, discuss the construction of possible correlation variables, and 

outline experimental challenges. We also briefly survey the Monte-Carlo tools which are useful 
\ in studying such correlations. 
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Most models of TeV scale new physics are motivated by solving the hierarchy problem. Therefore, 
the most crucial ingredient of all of them is the mechanism of cancelling the quadratically divergent 
correction to the Higgs mass within the Standard Model. In order to achieve this, a set of new 
physics particles with the same or similar gauge quantum numbers as the Standard Model particles 
are introduced, whose couplings to the Higgs are related to those of the Standard Model particles. 
This "partner" -like structure of new physics is very generic in large classes of new physics scenarios. 
Well-known examples include the set of superpartners in low energy supersymmetry [I] (for a recent 
review see Ref. [2]) , KK excitations in extra-dimensional models [3], as well as similar states in 
little Higgs models [4] (Ref. [5] provides a brief review). 

Due to the similarities in gauge quantum numbers, initial LHC signatures of new partners are 
very similar, as they can decay into the same set of observable final state particles. The mass spectra 
of different scenarios can be chosen to produce similar dominant kinematical features, such as the 
Pt distribution of the decay product. For example, a typical gluino decay chain in supersymmetry 
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is g — > qq + N2 followed by N2 — >• 11 + N\. A similar decay chain in universal extra-dimension 
models [3] with KK-gluon (gW), KK-W (W^ 1} ) and KK-photon (7W), ff M -► g<?IU 3 (1) followed by 

gives identical final states since both iVj. and 7 1 are neutral stable particles which 
escape detection. The mass spectra of both supersymmetry and UED can be adjusted in such a 
way that the px of the jets and leptons are quite similar. 

Some of the similarities in the LHC signature are actually the result of equivalences in low 
energy effective theory. For example, it is known that "theory space" motivated little Higgs models 
are equivalent to extra-dimensional models in which Higgs is a non-local mode in the bulk, via 
deconstruction [HUH [8]. Therefore, they can actually be described by the same set of low energy (~ 
TeV) degrees of freedom. An important feature of this class of models is that the partners typically 
have the same spin as their corresponding Standard Model particles. 

However, the difference between this set of models and low energy supersymmetry is physical 
and observable with a sufficiently precise measurement. In particular, the spin of superpartners 
differ from their Standard Model counter parts by half integers. Therefore, spin measurements are 
crucial to set these scenarios apart. 

The conventional way of measuring the spin of a new particle involves reconstruction of its rest 
frame using its decay products and studying the angular distribution about the polarization axis. 
For example, in process e + e~ — > Z — > the 1 + cos 2 9 distribution of the muon direction in 

the rest frame of the Z reveals its vector nature. However, in most new physics scenarios of interest 
such a strategy is complicated by the generic existence of undetectable massive particles. Motivated 
by electroweak precision constraints and the existence of Cold Dark Matter, many such scenarios 
incorporate some discrete symmetry which guarantees the existence of a lightest stable neutral 
particle. Well-known examples of such discrete symmetries include R-parity in supersymmetry, 
KK-parity of universal extra-dimension models (UED) [3], or similarly, T-parity in Little Higgs 
Models [9j [lOl QTJ [12] . The existence of such a neutral particle at the end of the decay chain results 
in large missing energy events in which new physics particles are produced. This fact helps to 
separate them from the Standard Model background. On the other hand, it also makes the spin 
measurement more complicated because it is generically impossible to reconstruct the momentum, 
and therefore the rest frame, of the decaying new physics particles. 

There are two different approaches to measuring spin. First, given the same gauge quantum 
numbers, particles with different spin usually have very different production rates, due to the 
difference between fermionic and bosonic couplings and the number of degrees of freedom. Such an 
approach could be useful, in particular initially, for colored particles due to their large (hence more 
measurable) production rates. However, a crucial ingredient in such a strategy is the measurement 
of the masses of particles produced, as rate can only provide definitive information once the mass 
is fixed. Such an effort is made more difficult owing to the existence of missing massive particles. 
There is also some residual model dependence since, for example, a couple of complex scalars can 
fake a dirac fermion. 

The second approach, is the direct measurement of spin through its effect on angular corre- 
lations in decay products. In the absence of a reconstructed rest frame, one is left to consider 
Lorentz invariant quantities which encode angular correlations. As we will see later in this review, 
spin correlations typically only exist in certain type of decays. Furthermore, new physics particles 
are frequently pair produced with independent decay chains containing similar particles. There- 
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fore, a valid spin correlation measurement requires the ability to identify a relatively pure sample of 
events where we can isolate certain decay chains and suppress combinatorics effectively. Therefore, 
except for very special cases, we expect this measurement will require large statistics. At the same 
time, as will be clear from our discussion, using the appropriate variables and correctly interpreting 
the measured angular distribution frequently requires at least a partial knowledge of the spectrum 
and the gauge quantum numbers. Obtaining information about the spectrum and the quantum 
numbers is likely to require a somewhat lower integrated luminosity than spin measurements do. 
Therefore, the order with which we uncover the properties of new particles is congruent to the or- 
der with which we must proceed in the first place to correctly establish these properties. Thus, we 
should clearly focus on mass scales, branching ratios and gauge quantum numbers first, once new 
physics is discovered at the LHC, while keeping an unbiased perspective towards the underlying 
model. More refined measurements, such as the ones described in this review, will enable us to tell 
the different models apart thereafter. Such measurements can be useful and even more powerful 
in a linear collider as was recently proposed in Ref. [13]. In this review we will concentrate on 
methods applicable to the LHC. 

In principle, the production of particles with different spins also leads to distinguishable an- 
gular distributions. This was investigated in the context of linear colliders in Ref. [14] . A similar 
measurement using the process pp — ► U* at the LHC has been studied in Ref. [15]. An analogues 
measurement in the production of colored states is more challenging. First, typically several differ- 
ent initial states and partial waves can contribute to the same production process. Therefore, it is 
difficult to extract spin information from the resulting angular distribution in a model independent 
way. Second, as commented above it is often difficult to reconstruct the direction of the original 
particles coming out of the production vertex. As a result, angular correlations are further washed 
out. 

In the rest of this review, we will survey both of these approaches with slightly heavier emphasis 
given to the angular correlation technique. For concreteness, we will compare supersymmetry with 
another generic scenario in which the partners, such as gluon partner g' , W-partner W, quark 
partner q' , etc., have the same spin as their corresponding Standard Model particles. As was 
pointed out above, this so called same-spin scenario effectively parameterizes almost all non-SUSY 
models which address the hierarchy problem. 

Spin measurement at the LHC is still a relatively new field where only first steps towards a 
comprehensive study have been taken. We will briefly summarize these developments in this review. 
We will focus here on the theoretical foundations and considerations relevant for the construction of 
observables. The potential for measuring spin in many new decay channels remains to be studied. 
Important effects, such as Standard Model background and large combinatorics, deserve careful 
further consideration. We outline these issues in connection to particular channels below. 



2 Rate and Mass Measurement 

The total cross section might serve as an initial hint to the spin of the new particles discovered 
[16] . Due to differences in the number of degrees of freedom, the couplings, spin-statistics and 
angular momentum conservation, particles with different spin have very different cross-sections. 
Two examples of such differences are shown in Fig. [TJ Therefore, if we could measure the mass of 
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Figure 1: Examples of the production rates of new physics particles. Top panel: the rate of Dirac 
fermion top partner vs right-handed stop (complex scalar) [T2]. Bottom panel: the rate of KK-gluon 
vs gluino (Majorana fermion) [16J. 



the new physics particles, it will provide us with a good measurement of the spin. Notice that many 
factors, such as the uncertainties in Parton Distribution Functions, total integrated luminosity, as 
well as the detector and cut efficiencies have to be taken into account to extract production rates at 
the LHC. Such difficulties could in principle be addressed with more detailed study and probably 
higher statistics. A sizable error bar on the production rate, say a factor of two, will not seriously 
affect these results since the rate in different scenarios are very different (see Fig [1]) due to the 
quantized nature of the spin. 

Certain model dependence is inevitable in this approach. For example, a fermion can be faked 
by two closely degenerate scalars. A potentially more serious problem is the model dependence 
in extracting production rates from observations. What we actually measure is the production 
rate multiplied by a some branching ratio into the particular set of final states we observe. Such 
branching ratios could easily vary by one order of magnitude over the parameter space of a certain 
scenario. Therefore, a straightforward extraction of rate information can only be achieved in the 
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cases where the branching ratio is simple and largely model independent. This could be possible 
since there are many examples in which a single channel dominates a one step decay process. It 
is also true when the coupling is simple and known, such as the decay of some colored particle. 
However, in more complicated cases and longer decay chains, more information is necessary to 
extract information about the production rate. 

Moreover, such a determination is only possible if we could measure the masses of the parti- 
cle using kinematical information. As demonstrated in p21[T7], typical "transverse" kinematical 
observables are mostly only sensitive to the mass difference between the decaying particle and the 
neutral particle escaping the detector. If this is all the information one can extract, there are 
clearly ambiguities in interpreting data in terms of an underlying model. For example, suppose 
supersymmetry is the correct model for TeV scale new physics, and we have observed gluino pro- 
duction at the LHC, pp — > gg, followed by g — > qq+ LSP. By assuming supersymmetry, we can find 
the gluino mass M g and LSP mass Mlsp which give rise to the measured rate and mass difference 
Am = M g — Mlsp- However, there is also, for example, an UED scenario in which KK-gluon with 
higher mass have the same production rate. It follows a similar decay chain g 1 — > qq+ LKP. We 
can adjust the mass of the LKP so that M g i — Mlkp ~ Am at which point this model results in the 
same experimental observables at the LHC. Notice that in order for this degeneracies to exist, the 
only requirement is that the observables are only functions of Am, even if the functions are different 
for different scenarios. Of course, we expect detailed measurement at the LHC to yield more than 
just the mass difference. With some assumptions regarding the underlying model there are more 
subtle kinematical observables which, in combination with the rate information, could determine 
the spin [17] . To what extent this could be generalized to a broader classes of new physics particles 
deserves further study. 

Recently, several methods have been proposed to measure the mass of new physics particles 
directly [13 [13 [201 EH [221 [231 E^- A detailed discussion of these methods is beyond the scope of 
this review. Combined with a rate measurement, such determination of mass will certainly provide 
an important measurement of the spin within the context of a certain (possibly well motivated) 
model. We expect such information to be complementary to the angular correlations we are about 
to discuss next. A detailed comparison between them with the effect of fully realistic Standard 
Model background and detector effects included has yet to be carried out. 



3 Angular Correlations in a General Decay Topology 

The decay topology we will be interested in is shown in Fig[2j A heavy partner, X, decays into a 
SM particle with momentum p\ and another partner Y . Y subsequently decays into a second SM 
particle with momentum p2 and a third heavy partner, Z which may or may not be stable. For most 
of the discussion we will assume that Y is produced on-shell, but it is important to keep in mind 
that this same topology may still apply even for 3-body decays where the intermediate particle is 
off-shell (we will explore this possibility in Section [7]). If Y is indeed on-shell, it is straightforward 
to show that the differential decay width is a polynomial in t\2 = (pi + P2) 2 of degree 2s, with s 
being the spin of Y (for details see |25| ) , 
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We have chosen to use the variable t±2 = mf 2 ) rather than invariant mass m\2- While they are 
equivalent in principle, distribution of t\2 is more convenient to interpret since it is linear in cos 6\2- 

The observed distribution often deviates from this simple relation because of both experimental 
and theoretical reasons. First, if the intermediate particle Y is not polarized, there is no hope of 
seeing any angular correlations, i.e. the polynomial is simply a constant. Second, if it is a fermion 
then polarizing it is not enough and we need its decay to distinguish handedness. This leads us to 
some simple necessary conditions for the existence of angular correlations [25J, 

1. If the intermediate particle Y is a vector-boson then My <C Mx so that it is longitudinally 
dominated and hence polarized. 

2. If the intermediate particle is a Dirac fermion then both vertices in the topology of Fig. [2] 
must be at least partially chiraQ- 

3. If the intermediate particle is a Major ana fermion then one must be able to determine the 
charge of the outgoing particles, p\ and p2 (even if only statistically). 



Pi 




Figure 2: The topology for a generic decay where spin information may be found. A heavy particle 
X decays into a SM particle with momentum pi and a heavy partner Y of spin s. Y subsequently 
decays into another SM particle with momentum p2 and an invisible particle Z. If Y is on-shell 
then the differential decay width for this process is a polynomial of order 2s in the relativistically 
invariant variable t\2 = (pi +P2) 2 - 

Based on these simple rules it is straightforward to list the different decay channels that may 
contain spin effects in the MSSM and other scenarios. In what follows, we divide these channels 
into those which probe the spin of the Electroweak gauge-bosons and Higgs sector partners and 
those which probe the spin of the matter sector partners. In certain cases additional information of 
the spin correlations can be used to determine the spin of the external particles, X and Z in Fig. 
[2j It can lead to a full spin determination of every new physics particle involved in a decay chain. 

1 We will see below that there is possibly an additional subtle kinematical suppression if the fermion consequently 
decays into a heavy vector boson, but this is non-generic. 
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4 Mis-pairing and Background 



Common to all attempts of establishing angular correlations is the problem of mis-pairing. The 
generic process shown in Fig. [2] is usually a part of a longer decay chain involving more final state 
particles. Moreover, this decay chain itself is only one of two branches in pair production. For 
example, the full event may look like Fig. [3] , which is a process commonly discussed in connection 
with spin determination and elaborated upon in Section ([5]). If we are interested in the spin of 
N2 we would like to correlate the adjacent outgoing jet and lepton, e.g. q\ and l near . However, 
in general there are no sharp kinematical differences between the near and far leptons along the 
same decay chain. Similarly, it is not easy to tell which is the correct jet to pair with the lepton. 
Therefore, an experimental histogram of the invariant mass variable t\2 is unlikely to follow the 
simple polynomial presented of Eq. (pQ). Instead we will get some shape which is composed of the 
polynomial structure coming from the correct pairing together with a more featureless distribution 
associated with the wrong pairing. This is in principle not a problem since one can fit to an 
underlying model by simulation, however, it does lead to a reduction in the relative statistical 
significance of the correct pairing. This problem is shared by all the channels we discuss below and 
we will refrain from the refrain where it is self-evident. 



Figure 3: A typical SUSY decay chain where the identification of the right pairing required for spin 
determination could be difficult. 

We will not discuss the problems associated with Standard Model background in this review 
and assume that it was reduced with a strong set of cuts on missing energy, number of leptons and 
etc. However, even if all the events contain new physics we may still have events with different 
underlying topology contributing the same visible final states. We refer to such ambiguities as same 
model background. Before investigating spin effects we must first ascertain that the event topology 
is well understood. Otherwise the resulting correlations will carry little significance. 



5 Spin Determination of Electroweak Gauge-Boson Partners 

In this section we cover potential channels which may reveal the spin of heavy partners of the EW 
gauge-bosons and scalars (in SUSY those would be the winos, binos and higgsinos). Considering 
the hierarchy problem it is necessary to have the gauge-boson partners at the EW scale since 
the gauge-boson loop contribution to the corrections to the Higgs mass are the largest after 
the top quark loops. We begin with the least model dependent channel and follow with channels 
which involve more assumptions about the spectrum. 
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5.1 Charged boson partner's spin - Set-W ± correlations 



The first channel we discuss is the minimal one in that it only assumes the presence of quark part- 
ners, charged EW bosons' partners and a stable lightest particle. This channel was first presented 
and investigated in Ref. [25J and later also in Ref. [26]. The relevant topology is shown in Fig. [J] 
for both SUSY and same-spin scenarios. 

The differential decay width for the two processes is a polynomial in t q w = (p q + pw) 2 , where 
p q is the momentum of the outgoing quark and pw is that of the outgoing W gauge-boson. It is of 
degree 1 (2) for the SUSY (same-spin) scenario and is given by, 

d^SUSY n ( , ,2 I i2\. . n dT same _ spin . 2 . . 

— = Cl{\aL\ - \CLR\ )t q w + 77 = <^2 t qW + ^\tqW + <^Q (2) 

atqW MqW 

where cll,r are the left and right coupling of the — C + — W + vertex and the Cj's are dimensionful 
factors which depend on the masses of the different particles in the decay and the quark partner's 
coupling. The leading coefficients are given by, 

_ g j M,(Ml-2Ml + Ml / 

° 2 = 256tt 2 M^i My M W i T w' ^ 

where (Xw) is the decay- width of the chargino (W) and we assumed that the quark partner 
coupling is dominated by its weak interaction. The expression for dT 'susy I 'dt q w in Eq.([2]) makes it 
evident that in the SUSY case angular correlations vanish in the non-chiral limit where \ai\ = \clr\- 




Figure 4: Topology for the decay of an up-type quark partner into a down-type quark, a W + 
gauge-boson and an invisible stable particle. On the left is the SUSY decay chain and on the right 
is the corresponding decay chain in a theory with same-spin partners. 



There are several experimental difficulties associated with this channel. First, if the W decays 
leptonically then one must correlate the resulting lepton with the outgoing jet. Since it is impossible 
to fully reconstruct the gauge-boson the distribution is no longer governed by the simple formula 
in Eq.([2|) (if the gauge-boson is strongly boosted in the lab frame one may be able to extract its 
four- vector at least approximately). Nonetheless, preliminary studies seem to indicate that the 
distinction between SUSY and same-spin scenarios is still very pronounced [25]. 
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Second, if the decays hadronically then it is important to establish how well can the 
be reconstructed and its momentum correlated with that of the outgoing jet. This will require a 
more realistic study and simulations. 

Fig-© contrasts the Jet-W^ correlation in SUSY against same-spin scenarios with quark 
partners at 1 TeV, vector boson partners at 500 GeV and an LSP at 100 GeV. This simulation 
was done including background from the same model (no SM background, though), i.e. events with 
identical final states, but different underlaying topology and without any knowledge of the correct 
pairing. In every event one tries to reconstruct the from two of the 4 jets and then form the 
invariant mass t q w with the two remaining jets. If more than one pairing reconstructs W it is 
regarded as failure and the event is discarded. The cuts involve > 200 GeV and \r]\ < 4.0. Jets 
were defined with a cone size of AR = 0.4. The linear behavior vs. the the quadratic behavior is 
still visible on top of the background coming from the wrong pairing and contamination from other 
event topologies. 



1.6 



1.5 



1.3 



Same-Spin 

SUSY 



Figure 5: Considering only signal events with 2-jets and a hadronic W this is a histogram of 
the invariant mass t q w (see text for details concerning the reconstruction) for the two prototype 
models: SUSY (solid-black), Same-Spin (dashed-red). The histogram is normalized to unit area 
and the normalized statistical error is approximately y/N ~ 0.04. The quadratic behavior of the 
same-spin scenario is clearly distinguished from the linear behavior associated with SUSY. 



This example also serves as a remainder that any observable suggested as a possible determinant 
of spin better result in robust differences in distributions' shapes. Any discriminator that relies on 
small discrepancies is unlikely to survive when the experimental limitations are properly taken into 
account. 
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5.2 Charged boson partner's spin - Jet-Z° correlations 

If the spectrum contains two light charged EW bosons' partners (e.g. both Higgsino and Wino are 
light) then the topology shown in Fig. [6] is possible. With sizable branching ratio, this is a very 
promising decay channel because the four-momentum of the Z° gauge-boson can be unambiguously 
reconstructed. The differential decay width follows that of Eq. ([2]) and the prospects of determining 
spin using this channel are now under investigation [27]. In addition, if the underlying model is 
indeed SUSY, the slope in Eq. (J2J) is sensitive to the value of tan/3. This is easy to understand 
if one recalls that for an intermediate Dirac fermion both interaction vertices need to be at least 
partially chiral. The q — C2 — q vertex is certainly chiral if the intermediate chargino is mostly a 
wino. The C2 — C\ — Z, however, is only chiral if tan/3 7^ 1. In fact it is straightforward to show 
that the slope of the distribution is directly related to tan/3 and vanishes as tan/3 — * 1. Since, this 
measurement is in principle very clean, it may offer extra information in addition to the spin of the 
chargino state. 




Figure 6: Topology for the decay of an up-type quark partner into an up-type quark, a Z gauge- 
boson and a charged partner. On the left is the SUSY decay chain and on the right is the corre- 
sponding decay chain in a theory with same-spin partners. 



5.3 Neutral boson partner's spin 

In an early and interesting paper [28], Barr has pointed out that the inherent charge asymmetry 
of a proton-proton collider may lead to certain angular correlations which carry spin information. 
This channel is known to be useful for mass determination and its sensitivity to spin was pointed 
out earlier in Ref. [29] and references therein. The decay topology for the SUSY case is shown in 

Fig. m 

In the case of SUSY, the intermmediate particle is a Majorana fermion and therefore, as pointed 
above, the charge of the outgoing quark and the near lepton must be determined in order to observe 
spin effects. The quark's charge may be obtained by either utilizing the charge asymmetry of the 
proton as pointed out by Barr, or by considering decays involving heavy flavor quarks such as 
bottom-quarks where a statistical determination of the charge is possible [30]. The lepton's charge 
can be determined unambigiously, but experimentally, the far lepton cannot be distinguished from 
the near lepton and one must average over both contributions. Barr found that it is possible to 
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define an asymmetry parameter which is sensitive to the spin of the intermediate neutralino N2, 

A _ ( dF _ dT \ I ( dT , dT 

\dt q l+ dtqi-J/ \dt q i+ dt q i- 

(The original asymmetry defined by Barr was with respect to the invariant mass rather than t q i±, 
but these two definition are isomorphic). As in the previous cases, this asymmetry has the opposite 
sign if the initial state is an anti-squark. Therefore, this observable relies on the LHC being a proton- 
proton collider and is sensitive to the quark - anti-quark asymmetries in the parton distribution 
functions. 

This decay chain was further investigated by [TH EU (32j [33] and contrasted against the cor- 
responding decay chain in the scenario with Universal Extra Dimensions (UED) j3J It is 
worthwhile to point out that when the intermmediate particle is a vector-boson (as in UED), no 
charge information is necessary in order to observe spin effects. However, correlations are washed 
out when the spectrum becomes degenerate. 





Figure 7: Topology for the cascade decay of a quark into two leptons and neutral stable partner. 
The correlation of the quark with the near lepton reveals the spin of the intermediate neutralino 
N\. However, experimentally it is hard to distinguish between the far and near lepton and one 
must average over both. The resulting distribution is asymmetric in the correlation of the quark 
with a positive versus negative leptons. 



We remark that this channel demands the existence of light lepton partners in the spectrum. 
Leptonic partners at the EW scale are certainly not required in order to solve the hierarchy problem. 
For example, they are not present in Little Higgs models. They exist in Supersymmetric theories 
and RGE effects tend to cause them to be lighter than the quark partners, but they may not 
necessarily be lighter than the neutralinos. Nonetheless, if they are light and the decay chain of 
Fig. (|7|) is allowed, it will teach us an enormous amount regarding the underlying model. 

Some of the experimental difficulties associated with this channel were pointed out above al- 
ready (mis-pairing and same model background) . In addition, the necessity for a charge asymmetry 
means that gluino pair production events will tend to contaminate the signal and reduce the sen- 
sitivity due to the absence of charge asymmetry to begin with. This is also true for squark pair 
production if dominated by gluon-gluon initial state. In connection with this difficulty, Ref. [30] 
had the interesting proposal of using the distinction between b and b. It may even lead to a mea- 
surement of the gluino spin with enough luminosity. However, so far this possibility has only been 
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investigated for the parameter point SPSla which is particularly friendly for hadron collider studies. 
It remains uncertain whether these conclusions can be extended to more generic spectra. 



6 Spin Determination of Standard Model Matter Partners 



In SUSY, the matter partners are scalars and therefore carry no spin information. In scenarios 
with same spin partners, angular correlations in decays involving the intermediate heavy fermions 
are certainly possible and the conditions for such effects were elucidated in Ref. [35J. It was also 
noted that, in combination with the measurements discussed above, one may obtain information 
on the spin of the initial and final partners in the cascade. We now briefly review these findings. 

The general decay topology is shown in Fig|8] where the SUSY decay of a gluino is shown 

alongside the corresponding decay in the same-spin scenario. The differential decay width is in 
general given by, 

d^SUSY dT same — S pi n . f 

~n— = c ° dt ff - = c ^f + c ° (6) 

where £// = (p/ + Pf) 2 and the coefficients C[ and Cq are presented in Table [2j 




Figure 8: Topology for a decay involving an intermediate matter partner. On the left is the SUSY 
decay of a gluino into two SM fermions and a neutralino. Since the intermediate particle is a scalar 
there are no angular correlations between the two outgoing fermions. On the right is the decay in 
a scenario with same spin partners. A heavy gluon-like particle decays into two SM fermions and 
a heavy neutral gauge-boson. The outgoing fermions do exhibit angular correlations owing to the 
fermionic nature of the intermediate heavy matter partner. 



As mentioned above, when the intermediate partner is a Dirac fermion the necessary conditions 
for the existence of angular correlations (i.e. for the coefficient C[ in Eq.([6]) to be non-zero) are that 
both interaction vertices in FigJS] are at least partially chiral. As we will now show, it is possible 
to satisfy this condition with very few assumptions about the same-spin model. 

First, we will assume that the members of this new sector have the same quantum numbers as 
the known quarks and leptons (this assumption is not necessary, but will simplify the discussion). 
As throughout this review, we will assume that some Z2 parity symmetry is present. In what 
follows we use the quark sector to demonstrate our results, but similar conclusions hold for the 
lepton sector as well. 



12 



SM 


Heavy partners 


SU{2) x U Y (l) 


Ql 

UR 

d-R 


Q'l = 


it) <* = 

U'r 
D'r 


V d' R J 
U' L 
D'l 




(1,1) 
(1,-1) 



Table 1: The matter content for the heavy partners in a generic same-spin scenario. 



We denote the heavy fermionic modes by Q' L R , U' L R and D' L R as shown in Table CO Here, 
Q' = (u',d') is an SU(2) doublet while U' and D' are singlets; note that Q' L , U' R and D' R are 
partners to SM fermions while Q' R , U' L and D' L have no Standard Model counterparts. 

At this stage, the coupling to the SM matter is via heavy bosons (</, W and etc.) because 
of the Z2 parity. Since the SM matter sector is chiral it forces the new interaction to be chiral as 
well. The coupling to the heavy gluon g' , for example, is schematically, 

Hint = Q'l^'il + U' R /u R + D' R ^'d R + h.c. (7) 

where qL is the SM electroweak doublet and ur and d R are the singlets. 

However, for the new fermions to be parametrically heavy they must have Dirac masses in 
addition to the usual Yukawa couplings to the Higgs. Therefore, after EW symmetry breaking, 
their mass matrix is given by (in here we ignore the flavor structure to simplify the discussion), 

Cmass,up ={W L U' L ) ( ^ ^ ) ( ^ ) • ( 8 ) 

where v = 246 GeV, A is a Yukawa coupling of the new states, and (p is some phase which in general 
cannot be rotated away. A similar mass matrix holds for the down sector. After diagonalizing the 
mass matrix the mixing angle is given by, 

2Xv JM% + 2M Q Mu cos ip + M% 

tan 6 = 2 ( 9 ) 

(m^ - A$) + J (Ml - + 4X 2 v 2 (m^ + 2M Q M V cos <p + M^j 

Below, we concentrate on two limiting cases of this relation which result in very different angular 
correlations. 



6.1 Non-degenerate Spectrum 

When the diagonal terms in the mass matrix are sufficiently separated, Mq — Mjj 3> Xv , the mixing 
between the states is small and the phase plays no role, 

tan9 ~XM~ ^ 1 (10) 

where AM = Mq — M\j. When the mixing is small the interactions with the SM in Eq.([7]) remain 
chiral after rotation into the mass eigenstate basis. This leads to an important conclusion: any 
model with heavy fermionic partners of the SM matter sector protected by some Z2 parity (KK- 
parity, T-parity etc.), with a non-degenerate spectrum, exhibits angular correlations in decays. 
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6.2 Degenerate Spectrum 



We now turn to examine a degenerate spectrum, so that Mq — Mjj <C Xv. Such a spectrum can 
result from any symmetry that relates left and right mass parameters. In this case the phase tp 
carries significance. When ip = the mixing is large, 

i AM 

tan6>= > 1 — (11) 

AM + a/ AM 2 + X 2 v 2 Xv 

The coupling of the mass eigenstates to the SM is no longer chiral, 

Ant = U'^'u + U'^'u + O (^f) + h.c. (12) 



and similarly for the other gauge couplings. Therefore, since the interactions with the SM, Eq. ()12p . 
are no longer chiral we expect angular correlations in decays to vanish (up to corrections of order 
0{AM/Xv)). 

In the UED model, not only are the masses degenerate by construction, but also the phase is 
fixed by 5-d Lorentz invariance ip = n as shown in Ref. [35]. In this case, unless A is unnaturally 
large, the mixing is always small, 

Xv Xv , . 

tan# = -> -= < 1 (13) 

M + \JM 2 + v 2 X 2 2M 

where M = (Mq + Mjj)/2. Therefore, as can be seen from Eq.([7]), the interactions of the SM with 
U[ and U!} remain chiral and angular correlations should be present in decays. 

This observation leads to the following conclusion in the case of degenerate spectrum: if angular 
correlations are not present, or are strongly suppressed, UED model can be ruled out. At the same 
time, it is still possible that matter sector partners are fermions. 



6.3 Slope Information 

If the vertices are chiral and angular correlations are therefore present then one may extract more 
information out of such measurements than just the spin of the intermediate matter partner. In 
Table [2] the other possible spin assignments for the external particles in a decay involving an 
intermediate Dirac fermion are showro All these distributions have an edge at, 

(m 2 -m 2 )(m 2 -m 2 ) 

However, depending on the spin of the external particles the slope is different and the behavior near 
the edge is modified. The slopes presented in the table assume no significant Left-Right mixing, 
such as UED or the non-degenerate case discussed above. A determination of the slope together 

2 The reader will undoubtedly notice that the slope vanishes as Mq —> 2My when 7' is a vector-boson. This subtle 
kinematic suppression arises because the outgoing vector-boson 7' is unpolarized in this limit. 
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with a measurement of the ratio Mg/Myi (from kinematical edges and cross-sections) can identify 
the spin of the external particle up to a two fold ambiguity. 

For example, if one measures a positive slope and Mg/My > 2 then one can conclude that 
the gluon partner is a vector-boson, but whether 7' is a scalar or a vector-boson is left unresolved. 
On the other hand, a positive slope together with Mq/M^, > 2 leads one to conclude that 7' is a 
scalar, but the spin of the gluon partner can be either zero or unity. 



Scenario 



Slope C[ 



Intercept C 




2M} 



M) 



M 2 r 



2M z gl 



r2 



M, 



2M 7 2 , 



Ml 



1 



{M? 



Q 



4M*M 2 g ,) t f} 



(edge) 



Ml t 



(edge) 
ff 



mK t 



(edge) 



Q "ff 



.(edge) 
ff 



Table 2: If angular correlations exist between the outgoing / — / or dilepton pair, then the sign 
of the slope of the distribution (whether C[ > or C[ < 0) may reveal the spin of the external 
particles as well as the intermediate one. In the first row we consider a scenario where the external 
particles are both vector-bosons (VFV). In the second row the incoming partner is a scalar whereas 
the outgoing partner is a vector-boson (SFV) and so forth. 



6.4 Long cascade decays and total spin determination 

The cascade decay in FigllOl is an example for a decay chain where the spin of all the partners 
may be determined unambiguously. Suppose one measures the slope of the bb pair to be negative 
with M|,/M§, < 2 and that of the dilepton pair, £~£ + , to be negative with M^/M 2 , < 2 as well. 
Then, either all three partners, g' , Z' and 7', are vector-bosons, or all three are scalars. Hence, 
with a single spin measurement of the Z 1 ', such as described in [2S\ UM I5T| 152"! I25| . one can lift 
this two- fold ambiguity and determine the spin of all the particles in the event. This optimistic 
picture is probably challenging in practice, however, it sharply illustrates the utility of the slope at 
distinguishing between different possible spin assignments. 

Measuring the spin in this way does not require reconstruction or a charge asymmetry. The 
most urgent challenge is to identify the correct pairing. Realistic simulations and experimental 
input will go a long way to shed some light on the issue. Also, the ability to tag the charge of the 
outgoing particles in the decays discussed above could be of use. A study which will quantitatively 
clarify how does the signal deteriorate when such information is only partially available (b tagging) 
or altogether missing (jets) is in need. This will be helpful for assessing the prospects for spin 
determination at the LHC if matter partners are light enough to be present in the spectrum. 
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Figure 9: The resulting differential decay-width distributions corresponding to the various possi- 
bilities discussed in Table [2j The mass scales are in units of M g >. The left pane corresponds to 
Mq/M-ji = 2 and the right one has Mq/M^i = 1.2. Notice the change in the slope for the VFV 
and SFV models as this mass ratio is changed. 




Figure 10: A long cascade decay in scenarios with same-spin partners. As explained in the text, a 
knowledge of the angular correlations between the bb pair and dilepton £~£ + may contain enough 
information to determine the spin of all the particles in the event, including g' , Z' and 7'. 



7 Off-shell decays 



As was discussed in Ref. [25J and elaborated upon in Ref. [36], if a decay proceeds through an off- 
shell fermion (consider FigJS] again, but with /' being off-shell), angular correlations are guaranteed 
because such a fermion is inevitably polarized. This can be understood by considering the fermionic 
propagator for a Dirac fermion, 

Mp) = htp (15> 

For M S> ft, the propagator favors one helicity structure over the other or in other words the 
intermediate particle is polarized. In SUSY, the intermediate particle is a scalar and no correlations 
are expected. The decay's kinematics are then governed by 3-body phase-space and the resulting 
differential decay-width is shown in Figllll This is contrasted with the same-spin scenario where the 
polarized intermediate fermion affects the distribution and causes it to deviate from the phase-space 
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prediction. 




Figure 11: Theoretical curves for the differential decay-width in decays of the type shown in Figj2]as 
a function of ti2 = (pi +P2) 2 for an off-shell intermediate state. The graphs are normalized to unit 
area. We contrast SUSY against two type of same-spin models. In SUSY, the intermediate particle 
is a scalar whereas the decaying particle and the LSP are Majorana fermions i.e. the gluino and 
neutralino. In VFV (SFS) scenario, the intermediate particle is a Dirac fermion and the external 
partners are vector-bosons (scalars). 

A simple example of such a decay is realized in scenarios where the matter partners are heavy, 
but the gluon partner can be produced directly (for an extreme realization of such a setup see 
Ref. [37]). It will then decay into neutral matter through an off-shell matter partner as in Fig. [8] 
(though, in the case of an intermediate off-shell state, one must be careful to include the crossed 
diagrams as well). 

On top of the experimental difficulties mentioned above, there is one additional difficulty that 
must be recognized. A brief look at the distributions in Figllll reveals that the real discrepancy 
between the different scenarios is concentrated in the low t±2 region. This corresponds to the 
outgoing quarks being close by in real space. It is therefore susceptible to AR and isolation cuts 
and it is not even clear that such quarks will be tagged as two different jets by the jet algorithms. 
These effects can carve out the entire low ti2 region and reduce the sensitivity of this observable 
to spin information. These effects have been included in the study of Ref. [36], but since no jet 
algorithm was used, it would be important to reanalyze this channel with a more realistic simulation. 
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7.1 Simulation tools to study spin correlations 



The amplitude for a single process is represented at tree level by several topologically distinct 
Feynman diagrams. The denominator for the amplitude corresponding to each diagram comes 
in the form of a product of simple poles. The numerator on the other hand takes into account 
the helicity structure of the diagram and can be written as the sum of several terms. To obtain 
an accurate simulation of the process one must square the full amplitude and include all the 
interference terms. Simulators such as CompHEP [38] and MadGraph [39J achieve exactly that 
albeit through slightly different routes (CompHEP offers in addition an analytic expression for the 
squared amplitude). 

However, when interested in events with multiple final states (more than 4 or 5), it becomes too 
computationally demanding to proceed without an approximation. If the intermediate states within 
each diagram can go on-shell they will do so since such configurations receive strong support from 
the pole structure of the denominator. In this case, the interference terms between denominators 
with different poles are suppressed with respect to the leading singularities. It is then possible 
to approximate the full amplitude by a reduced one in which the intermediate states are set on- 
shell and any interference is neglected. This is the well-known narrow-width approximation and is 
utilized by the most widely used generators: PYTHIA [ID] and Herwig |41j . This approximation 
of course breaks down near a threshold or when the width becomes large compared with the mass 
difference (in which case the notion of an isolated pole identified as a particle is meaningless to 
begin with). It goes without saying that in the case of an off-shell decay such an approximation is 
altogether invalid. 

As we just saw that the narrow-width approximation ignores interference terms between topo- 
logically distinct Feynman diagrams. There may still be additional interference effects associated 
with the different polarizations of intermediate particles in individual Feynman diagrams. PYTHIA 
and Herwig differ in the way they treat these effects. In consequence of the narrow-width approx- 
imation the propagator of intermediate states is reduced to the sum over polarization and we can 
describe a generic matrix element involving one propagator schematically as, 



A 

where the sum is over the different intermediate polarizations A. The event-generator is interested 
in the squared amplitude which corresponds to the probability, 



In the last stage we defined the spin density matrix D\y which embodies all the information carried 
by the numerator. Now, the spin of the initial and final state is presumably known either because 
we must average or sum over it. The question is what is being assumed about the density matrix 
and whether it is to be calculated. PYTHIA simply ignores it all together setting it to the identity 
matrix properly normalized. Herwig uses a clever algorithm and computes it stochastically [42\ 143] . 
It therefore maintains full coherence as far as the spin is concerned. It was recently extended to 
include heavy on-shell vector-bosons (such as KK gluon states etc.) and the details are described 




(16) 





(17) 



in GS 



18 



As mentioned above the problem in using fully coherent event generator is that it is computa- 
tionally prohibitively expensive when the number of external particles exceeds a few. Events which 
result in many final states are unwieldy using this scheme. Recently, a new and useful tool called 
Bridge [44] was developed to deal with longer decay chains by making use of the narrow- width ap- 
proximation. Unlike PYTHIA, Bridge does not ignore the density matrix altogether, but rather try 
to approximate it by computing the diagonal elements, D\\, in a particular reference frame. This 
results in a fairly good approximation in certain cases when one of the elements is much larger than 
the others (such as in Top decay through an intermediate W-boson). However, it should not be 
expected to hold in general. In particular, since the density matrix is not Lorentz invariant, while 
the approximation may be valid in one reference frame it is not necessarily so in all frames. Which 
frame is optimal cannot be answered in full generality. It may prove useful to try and incorporate 
the algorithm for maintaining spin coherence used in Herwig into Bridge. 



8 Conclusion and Outlook 



We have reviewed methods of measuring the spin of new physics particles at the LHC. Such mea- 
surements are crucial in distinguishing qualitatively different classes of new physics scenarios, such 
as supersymmetry and Little Higgs or extra-dimensional models. 

For given masses, new physics particles with similar gauge quantum numbers but different 
spin will typically have very different production rates. Therefore, combining with independent 
measurements of their masses, their production cross-sections provide a quick determination of 
their spin, with minimal model assumptions. However, due to the existence of neutral stable 
massive particle at the end of the decay chain in many new physics scenarios, absolute mass 
measurement is challenging. Moreover, model dependent factors, such as decay branching ratios, 
have to be taken into account to extract total production rate. New method of utilizing the rate of 
information, such as combining production rates in several channels and using certain kinematical 
features which are sensitive to the overall mass scale, needs to be developed. Several new approaches 
of measuring absolute mass scales have been proposed recently. Their effectiveness in providing 
additional information in spin measurement need to be accessed. 

We focus on direct and model independent ways of measuring spin using angular correlation 
among decay products of the new physics particles. Due to the difficulties in fully reconstructing of 
the rest frame of the new physics particle from its decay products, we have to use Lorentz invariant 
combinations of momenta of the observable particles, which encode the spin information of the 
intermediate particles in the decay chain. We survey the potential of extracting spin information in 
a large class of different decay topologies. In each case, we present specific observables and discuss 
the conditions under which spin correlations are in principle observable. 

We considered two classes of new physics particles: partners of the Standard Model gauge- 
bosons and partners of the Standard Model matter fermions. In both cases, understanding the 
masses of particles and the chiral structure of their couplings is crucial in extracting the spin 
information. Details of exclusive decay chains must be well understood so that the appropriate 
variables can be used and interpreted correctly. As illustrated in the text these issues apply directly 
to the gauge-boson partners. In the case of fermionic Standard Model matter partners The couplings 
are always chiral when the spectrum is non-degenerate. In the more subtle degenerate case the 
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existence of spin effects depends on additional parameters. In particular we concluded that the 
UED model always leads to angular correlations as a result of five-dimensional Lorentz invariance. 
Moreover, the slope of the resulting distributions contains additional information regarding the 
spin of the external particles. Its determination requires knowledge of the masses of those particles. 
Building on these findings, we also propose a way of measuring in principle the spin of every new 
physics particles in the decay chain. In the case of an off-shell decay with an intermediate fermion, 
angular correlations are always present. This is an interesting direction which deserves further 
detailed investigation. 

Both the existence and the interpretation of the spin correlation at the LHC depend strongly on 
the parameters of the model under consideration. Therefore, we need to explore potentially useful 
channels and observables as much as possible. Viable strategies of extracting spin information at 
the LHC can only be finalized with our measurements of masses and couplings of the new physics 
particles. It is important to keep as much viable frameworks in mind as possible while we only 
have incomplete information about the new physics discovered at the LHC. 



References 

[1] Savas Dimopoulos and Howard Georgi. Softly broken supersymmetry and su(5). Nucl. Phys., 
B193:150, 1981. 

[2] D. J. H. Chung et al. The soft super symmetry-breaking lagrangian: Theory and applications. 
Phys. Rept, 407:1-203, 2005. 

[3] Thomas Appelquist, Hsin-Chia Cheng, and Bogdan A. Dobrescu. Bounds on universal extra 
dimensions. Phys. Rev., D64:035002, 2001. 

[4] Nima Arkani-Hamed, Andrew G. Cohen, and Howard Georgi. Electroweak symmetry breaking 
from dimensional deconstruction. Phys. Lett, B513:232-240, 2001. 

[5] Martin Schmaltz and David Tucker-Smith. Little higgs review. Ann. Rev. Nucl. Part. Sci., 
55:229-270, 2005. 

[6] Nima Arkani-Hamed, Andrew G. Cohen, and Howard Georgi. (de) constructing dimensions. 
Phys. Rev. Lett, 86:4757-4761, 2001. 

[7] Hsin-Chia Cheng, Christopher T. Hill, Stefan Pokorski, and Jing Wang. The standard model 
in the latticized bulk. Phys. Rev., D64:065007, 2001. 

[8] Hsin-Chia Cheng, Christopher T. Hill, and Jing Wang. Dynamical electroweak breaking and 
latticized extra dimensions. Phys. Rev., D64:095003, 2001. 

[9] Hsin-Chia Cheng and Ian Low. Tev symmetry and the little hierarchy problem. JHEP, 09:051, 
2003. 

[10] Hsin-Chia Cheng and Ian Low. Little hierarchy, little higgses, and a little symmetry. JHEP, 
08:061, 2004. 

[11] Ian Low. T parity and the littlest higgs. JHEP, 10:067, 2004. 



20 



[12] Hsin-Chia Cheng, Ian Low, and Lian-Tao Wang. Top partners in little higgs theories with 
t-parity. 2005. 

[13] Matthew R. Buckley, Hitoshi Murayama, William Klemm, and Vikram Rentala. Discriminat- 
ing spin through quantum interference. 2007. 

[14] Marco Battaglia, AseshKrishna Datta, Albert De Roeck, Kyoungchul Kong, and Konstantin T. 
Matchev. Contrasting supersymmetry and universal extra dimensions at the clic multi-tev e+ 
e- collider. JEEP, 07:033, 2005. 

[15] A. J. Barr. Measuring slepton spin at the lhc. JEEP, 02:042, 2006. 

[16] AseshKrishna Datta, Gordon L. Kane, and Manuel Toharia. Is it susy? 2005. 

[17] Patrick Meade and Matthew Reece. Top partners at the lhc: Spin and mass measurement. 
Phys. Rev., D74:015010, 2006. 

[18] Hsin-Chia Cheng, John F. Gunion, Zhenyu Han, Guido Marandella, and Bob McElrath. Mass 
determination in susy-like events with missing energy. JEEP, 12:076, 2007. 

[19] Won Sang Cho, Kiwoon Choi, Yeong Gyun Kim, and Chan Beom Park. Gluino stransverse 
mass. 2007. 

[20] Ben Gripaios. Transverse observables and mass determination at hadron colliders. 2007. 

[21] Alan J. Barr, Ben Gripaios, and Christopher G. Lester. Weighing wimps with kinks at colliders: 
Invisible particle mass measurements from endpoints. 2007. 

[22] Won Sang Cho, Kiwoon Choi, Yeong Gyun Kim, and Chan Beom Park. Measuring superpar- 
ticle masses at hadron collider using the transverse mass kink. 2007. 

[23] Graham G. Ross and Mario Serna. Mass determination of new states at hadron colliders. 2007. 

[24] M. M. Nojiri, G. Polesello, and D. R. Tovey. A hybrid method for determining susy particle 
masses at the lhc with fully identified cascade decays. 2007. 

[25] Lian-Tao Wang and Itay Yavin. Spin measurements in cascade decays at the lhc. JEEP, 
04:032, 2007. 

[26] Jennifer M. Smillie. Spin correlations in decay chains involving w bosons. Eur. Phys. J., 
C5L933-943, 2007. 

[27] A. Hook, L. Wang, and I. Yavin. to be published. 

[28] A. J. Barr. Using lepton charge asymmetry to investigate the spin of supersymmetric particles 
at the lhc. Phys. Lett, B596:205-212, 2004. 

[29] Peter Richardson. Spin correlations in monte carlo simulations. JEEP, 11:029, 2001. 

[30] Alexandre Alves, Oscar Eboli, and Tilman Plehn. It's a gluino. Phys. Rev., D74:095010, 2006. 

[31] Jennifer M. Smillie and Bryan R. Webber. Distinguishing spins in supersymmetric and uni- 
versal extra dimension models at the large hadron collider. JEEP, 10:069, 2005. 



21 



[32] AseshKrishna Datta, Kyoungchul Kong, and Konstantin T. Matchev. Discrimination of su- 
persymmetry and universal extra dimensions at hadron colliders. Phys. Rev., D72:096006, 
2005. 

[33] Christiana Athanasiou, Christopher G. Lester, Jennifer M. Smillie, and Bryan R. Webber. 
Distinguishing spins in decay chains at the large hadron collider. JHEP, 08:055, 2006. 

[34] Hsin-Chia Cheng, Konstantin T. Matchev, and Martin Schmaltz. Bosonic supersymmetry? 
getting fooled at the lhc. Phys. Rev., D66:056006, 2002. 

[35] Can Kilic, Lian-Tao Wang, and Itay Yavin. On the existence of angular correlations in decays 
with heavy matter partners. JHEP, 05:052, 2007. 

[36] Csaba Csaki, Johannes Heinonen, and Maxim Perelstein. Testing gluino spin with three-body 
decays. JHEP, 10:107, 2007. 

[37] Nima Arkani-Hamed and Savas Dimopoulos. Supersymmetric unification without low energy 
supersymmetry and signatures for fine-tuning at the lhc. JHEP, 06:073, 2005. 

[38] E. Boos et al. Comphep 4.4: Automatic computations from lagrangians to events. Nucl. 
Instrum. Meth., A534:250-259, 2004. 

[39] Fabio Maltoni and Tim Stelzer. Madevent: Automatic event generation with madgraph. JHEP, 
02:027, 2003. 

[40] Torbjorn Sjostrand, Stephen Mrenna, and Peter Skands. Pythia 6.4 physics and manual. 
JHEP, 05:026, 2006. 

[41] G. Corcella et al. Herwig 6.5 release note. 2002. 

[42] John C. Collins. Spin correlations in monte carlo event generators. Nucl. Phys., B304:794, 
1988. 

[43] I. G. Knowles. Spin correlations in parton - parton scattering. Nucl. Phys., B310:571, 1988. 

[44] Patrick Meade and Matthew Reece. Bridge: Branching ratio inquiry / decay generated events. 
2007. 



22 



